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X. 


Introduction 


The Time/Data 100 is a small digital computer specialized in 
time- series analysis. 

Its algorithms include the Direct Fourier Transform, the Auto- 
Spectral Density and the Autocorrelation. 

This report will analyze the Power Spectral Averaging Algorithm 
[l]. This algorithm applies the Direct Fourier Transform Algorithm to 
an input record. Then it squares and adds the real and imagina / 
components of the Direct Fourier Transform and accumulates each successive 
output record. 

An input record is a 1001 wor< block, where each word has 8 bits. 

An output record also consists of 1001 words, but has 18 bits per word. 

The computer has an oscilloscope and a plotter for analog output 
and tapes for digital output but they can only display 8 bits of the 
18 bit output word. 

The way the Direct Fourier Transform Algorithm operates and the 
limited display range of an output word dictated a need for a scale 
factor so that the output data could be analyzed properly. 

A scale factor was derived by analyzing the process of computing 
the power spectrum of an input data record. Each part of the scale 
factor is examined and explained as to why it's needed. The last 
section will show some examples to verify the above scale factor. 
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II, Input Data Factor 

Fo is the allowable maximum positive value of the input data, 

Fo “ F/ (data normalization factor) 

where F « 127 * 2 7 ~1 is the maximum number that the Time/Data 100 can 
receive in an 8 bit input word. 

The data normalization factor is a number that is multiplied to 
the data before it is used as input to the Time/Data 100. It is used 
to scale the data so that all of it will be within the + 127 range. 

Note that the Time/Data 100 treats the total range as + 1 therefore 
the data is automatically divided by 127. The factor F arises from 
this fact. For example, in our case the data values were so small 
that we could multiply them by 10 to achieve greater accuracy. 

Therefore Fo = F • (1/10) =» 12,7 

Fo is squared because the operation was completed before execution 
of the Spectral Averaging Algorithm. Therefore FcP must be multiplied 
to the output data. 
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XXX, Input Scale Factor 

Nl Is the input scale factor. 

The Time/Data 100 outputs an 18 bit word for each point in the 

Direct Fourier Transform Algorithm. Since the Spectral Averaging 

Algorithm only uses 7 bits plus the sign bit of the 17 bits plus the 

sign bit output word for its input, the operator must select a 7 bit 

range with the amplitude scaling switch that will include the most 

significant bits of the 18 bit output word. The selected number on 

the amplitude scaling switch corresponds to Nl. This, in effect, 

N1 

divides the original word by a factor of 2 ; Nl being the lowest 

order bit number of the 18 bit word that the 7 bits were selected 

from. For example, if the most significant bits were from 2 s through 

Z 11 , then by setting the amplitude scaling switch to 5 that bit range 

will be chosen for input to the Spectral Averaging Algorithm. So 

the original word is divided by 2 B . 

Since this is also done before the Spectral Averaging Algorithm, 

2N1 , 

the factor is squared. So, 2 must be multiplied. 




IV. Output Scale Factor 
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N2 is the output scale factor. 

The Spectral Averaging Algorithm also outputs a 17 bits plus the 
sign bit word, but the analog and digital output devices can only 
display 7 bits plus the sign bit of a word. So, a dynamic range must 
also be selected on the amplitude scaling switch to display the most 
significant bit. As in the input scale factor, this has the effect 
of dividing the output word by 2 N2 . So, 2 N2 must be multiplied to 
counteract the previous division. 

V. Data Shifts 

The Time/Data 100 uses an implementation of the Rapid Fourier 
Transform for the Direct Fourier Transform Algorithm. While trans- 
forming the data, some data folds are used. Of these data folds, two 
are sums and they produce 9 bit output words, 8 bits plus a sign bit. 
To conform to the multiplier requirements, which can only take 7 bits 
plus a sign bit, the data must be shifted down one bit. This is done 
twice, so the data is divided by 2 s . This is also done before the 
Spectral Averaging Algorithm, so 2* must be multiplied to the output 
data. 


VI. Averaging Factor 

N' is the number of output records that have been accumulated. 
The Spectral Averaging Algorithm only adds the output records. So, 
to obtain an averaged spectrum, the output data must be divided by 


the number of accumulated records. 
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VII, Algorithm Implementation Factor 


The final form of the scale factor had to be derived by analyzing 
the Time/Data 100’s method of computing the average value of the 
power spectrum and comparing it to the true value of the power spectrum. 

Given a time X, suppose that there exists a continuous function 
F(t 4- X) , *»<t< M , that satisfies the ergodic hypothesis and a square 
window B(t) such that 

B(t) - 1 for |t| & Tn/2 

B(t) •• 0 for |t| > Tn/2 

B(t) “ B (-t) 

for some interval of time Tn. 

Let h(t + X) * F(t + X) B(t) , -Tn £ t £ Tn, Partition the interval 

-Tn £ t £ Tn into 4N + 1 parts, where N is any integer such that 

N ^ 0 and let At be the length of one partition such that (4N -1- 1) At * 2Tn, 

Finally, let the discrete function g(k), -2N £ k £ 2N be the set of 

points of h(t) such that they are the midpoints of each partition. Figure ] 

describes the partitioned interval. Note that the distance between 

two points of g(k) is also At, which is the sampling rate. 


Figure 1, 
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(4N + 1) At = 2Tn 


This midpoint of each partition is a point of g(k). 


The Direct Fourier Transform Algorithm assumes that a 1001 point 
input record is appended by 500 zero points on each side. Hence N 
is 500 in this algorithm. 

So, the Time /Data 100 Computes the Direct Fourier Transform with 
this approximation? 

v ~12rr (jAf) (mAt) 

A(.f) ■ . f? (to) e (1) 

irJ 

m 


(Throughout this report the variables j, 1, and m vary from -2N to 2N 

< + X 2N 

and means 2, unless otherwise specified.) 

«2N 

where At is the sampling rate and Af “ 1/2 Tn * 1/ (4NF1) At 

~i2njm/(4NFl) 

A(j) * ^ g(m) e 
m 


The Power Spectral Density Algorithm then computes: 
P(j) - A(j) • A* (j), 


where * represents the complex conjugate. 


*(J> 


m 


t 


g(m) g(r) e 


i2rrj (m«r)/ (4NF1) 


This approximation, however, doesn't use the sampling rate, At, 
or the length of the input interval, Tn, in the computation of the power 
spectrum. 


Rewrite equation (1): 

A(jAf) - (1/At) 'l g(fflAt) e * i2rr(jAf) 

m 

Express the summation as an integral and replace the discrete 


function g(mAfc) by its continuous counterpart F(t+\) B(t) 
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A (£,,X> - (1/At) | F (t+X) B(t) a" 1217 *^ dt 

where the integrals have the limit* - » to ®; this applies to what 
follows. 

The Power for frequency £j is: 

P(£j»X) - (1/At) a J'JV (t+X) B(t) F(s+X) B(s)e i2nf j ( ' t " 8 ^dBdt 


Ensemble averaging is defined by: 


P r (£j) 


r* T 


* lim 1_ 

T -* » 2T _ T 


P(£j,X) dX 


There P’(fj) is the average value of P(fj,X) over all possible 
records. 


'■ < V 


T 

lim (1/2T)P (1/At) 2 IT F(t+X) B(t) F(s-fX) B(s)e 
T -» 03 { T 


i2nf . (t-s) 

^ dsdtdX 


Rearrange Terms : 


P'(£j) = (l/At) 2 Jj B(t) B(s)e 


i2rrf ^ (t~s) 


lim (1/2T)J F(s+X) F(t+X) dXdsdt 
T 00 -T 


Let T = t + X 

dT = dX and 8*X = (s-t) +T 
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P'(fj) «(1/At ) ? jJ B(t) B(s) e i2T7f j (t “* ) lim (1/2 T)|f[ (s-O+t] F[T]aTdsdc (2) 


Note that the last expression is the autocovariance of F(t). 

I.et CU-t) “ Hm (1/2T)JF [(a-t) + t] F(t) dT 

*1* — * CO ~T 

The autocovariance function is also defined by: [2] 

C (o) ■ | G (f) e* 217 ^ df W 

where 

G (f) * J C (a) e da 

is the power spectral density function. 

Substitute (3) into (2) 

P r (fj) * (1/At) 3 JJ* B Ct) B(s) e i2rrf j^ t 8 ^J* G(f) e dfdsdt 


P' Cf j) - (1/At? ) J G(f)dfJ’B(t)e i2TTt(f j‘ f) dtJ’ B(s)e" i2r,s(f j" f) ds 
P' (f ) - <1/At?) J* G(£) H(f ,-£) H* (fj-f) df 

J 

where H (f) * Jb (t) e ^^ rr ^ t dt 

and * represent the complex conjugate 

P' (fj) * (1/A t 8 ) J G(f) | H(fj-f) | s df 


Except for the scale factor, 1/At 3 , this is the convolution of 
the power spectral density with a spectral window. 

Compute the spectral window, 

H (fj-f) “ J B(t) e i 2 rr ( f j- f ^ 

,Tn/2 12n(f r f)t dt 
" -Tn/2 

_ / irr(f j -f)Tn _ e -iTT(f r f)Tn^ / i2lT ( fj .f) 


( 4 ) 
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*■ Tn(co» tt0 + i Bin tt0 - £cob rr0 - i ain rr0]) / I2n0 

where 0 * (fj-f) Tn 

H(f -f) - Tn sin (rr0)/rr9 - Tn dif 0 
j 

Rewrite equation (4) 

Ac® P'Cfj) - J* G(f) In® (dif 6)® df 
(Ac®/Tn®) P'(fj) “ J“ G(f) (dif 9)® df 


(5) 


The derivation, in effect, stops here. So the spectrum must 
be interpreted with a spectral window of (dif 0) s . The window is 
pictured in Figure 2. 

For the Time/Data 100, we have to replace P'(fj) with the machines 
estimate along with the previously defined factors. 

P'(fj) - (f<^ 2 2N1 + N2 + 4 /N') P td (fj) 

where P^ D (fj) is the machine's power spectral estimate. 

Rewriting equation (5) : 

[Fcf 3 2 2Nl + N2 +6 / N’ (4N 4-1) 3 ] VtjrQ (f j) * Jg (f) (dif 9) 3 df (6) 

where At 3 /Tn 3 ** 4/ (4N + l) 3 

This can be simplified if the power has a smoothly varying spectrum 
and G(f ) is relatively constant around f^. Then, we can approximate 
equation (6) by: 

(fo s 2 2N1 + 1,2 + 6 / N* (4N + 1)* P m (fj) s* G(f j) 6f 
where 6f ■ J (dif 0) 2 df ** 1/Tn 
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([fc? 2' £N1 + N2 + S At] / N' (4N + 1)) P TO (fj) *> G (f ) 

This estimate is, however, for a double sided spectrum. Therefore, 
to obtain a single sided spectral estimate, it must be multiplied 
by 2, namely 

(fo* 2 2N1 + 1,2 + 6 At/ N’ <4N + 1)) P to (fj) w G (£j). 
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VIII. Examples and Conclusion 

To test the scale factor, some examples are taken from the 
rubidium magnetometer experiment on the 0G0-5 satellite that measures 
the magnitude of the intensity of the magnetic field. 

Two methods were used to analyze the test data. First, by using 
the Time/Data 100 Method and second by using the Blackman-Tukey 
Method of Power Spectral Analysis. 

First, a number of records were accumulated on the Time/Data 100 
and the factors were noted. The output was then plotted on the output 
plotter. 

The same data set was then processed by computer with the Blackman- 
Tukey Method. To compare the two methods, the scale factor was computed 
and the output from the Blackman-Tukey Method was plotted according 
to the scale set by the computed scale factor. 

To test all portions of the scale factor, the data was sampled 
at three different rates; this change in At brought about a 
corresponding change in Nl, N2, and N 1 of the scale factor. 

First, the data was sampled at 1 point per 144 ms and graphs 1 and 
2 show the results of the Time/Data 100 method and the Blackman-Tukey 
method respectively. 

Graphs 3 and 4 show the two sampled at 1 point per 1.008 seconds 
and 5 and 6 show them sampled at 2.016 seconds. 

i 

Upon comparison, these results clearly show that the scale factor 
does produce the correct results. 
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